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Exact solutions are derived for the problem of a two-dimensional, infinitely anisotropic, linear-elastic medium 
containing a periodic lattice of voids. The matrix material possesses either one infinitely soft, or one infinitely 
hard loading direction, which induces localized (singular) field configurations. The effective elastic moduli are 
computed as functions of the porosity in each case. Their dilute expansions feature half-integer powers of the 
porosity, which can be correlated to the localized field patterns. Statistical characterizations of the fields, such as 
their first moments and their histograms are provided, with particular emphasis on the singularities of the latter 
The behavior of the system near the void close packing fraction is also investigated. The results of this work 
shed light on corresponding results for strongly nonlinear porous media, which have been obtained recently 
by means of the "second-order" homogenization method, and where the dilute estimates also exhibit fractional 
powers of the porosity. 

PACS numbers: A voids and inclusions; B anisotropic material; B constitutive behaviour; C Energy methods; localization. 



I. INTRODUCTION 
A. Context 

Most of the currently available homogenization methods for strongly non-linear composites make use an underlying linear 
homogenization estimate, whether aimed at computing dielectric (Willis, 1986; Zeng et al., 1988; Ponte Castaneda, DeBotton, 
Li, 1992), or elastic-plastic transport properties (Ponte Castaneda, 1991; Ponte Castaneda, 1996; Masson et al., 2000; Ponte 
Castaneda, 2002; Lahellec and Suquet, 2004). The best results are obtained using an anisotropic linear estimate, the anisotropy 
of which is consistently determined, by the means and variances of the fields in each phase of the composite (Pellegrini, 2001; 
Ponte Castaneda, 2001). The anisotropy of the underlying linear medium accounts for the privileged direction imposed by the 
driving field in the non-linear medium. 

Applying the "second-order" method (Ponte Castaiieda, 2002) to a random power-law material weakened by aligned cylin- 
drical voids, it has been found that in the so-called "dilute" limit (i.e., that of a vanishingly small volume fraction of porosity 
/), and in the limit of infinite exponent where the matrix becomes ideally plastic with a yield threshold, the leading correction 
to the yield stress in pure shear behaves as f^^"^. This result holds for both Hashin-Shtrikman and self-consistent estimates. 
The second-order method of Pellegrini (2001), suitably modified by replacing a Gaussian ansatz for the field distributions by a 
Heaviside ansatz (so as to cut-off high fields in a way compatible with threshold-type materials), provides in some cases similar 
predictions (Pellegrini and Ponte Castaiieda, 2001, unpublished). This /^/"^ dilute behavior is unusual in the context of effective- 
medium theories for random media, where the low-order terms in dilute expansions from the literature usually consist in integer 
powers of /. With the f^^^ correction, the derivative of the yield stress with respect to the porosity is infinite at / = 0, due to the 
exponent being lower than 1 . Thus, a vanishingly small volume fraction of voids induces a dramatic weakening of the porous 
medium. Ponte Castaiieda (1996, 2002) interprets this phenomenon as an indication of localizing behavior in a regime where 
shear bands pass through the pores, remarking that the limit analysis of Drucker (1966) produces an upper bound for the yield 
threshold with a dilute correction behaving as /^/^ in 2D. Numerical calculations and tests on perforated plates with periodically 
distributed holes (Francescato and Pastor, 1998) are consistent with these predictions in a plane stress situation. However, for 
a square network of circular holes, kinematic and static limit analyses result in the following analytical bounds on the effective 
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yield stress Y (Francescato et al., 2004): 



1 - 2(//7r)i/2 < Y/Y < (2/V3) [l - 2(/A)i/2 



(1) 



where Y is the yield threshold of the matrix. Numerical limit-analysis on a hollow disk model in plane strain leads to somewhat 
different conclusions, with a correction behaving roughly as /^/■^ for uniform strain boundary conditions, and as /^/^ for uniform 
stress boundary conditions (Pastor and Ponte Castaneda,2002). However, the calculations were not sufficiently accurate to be 
completely definitive. On the other hand, the exponent is quite certainly associated with localizing behavior, and its presence 
seems to be independent on whether the system is random, periodic, or comprises a unique void (although the actual values of 
the exponent may be dependent on the specific microstructure considered). 

Because the second-order theory with exponent 2/3 inherits its properties from an underlying anisotropic linear effective- 
medium theory, whose anisotropy is consistently determined by the field variances, a natural question concerns the "localizing" 
properties of the anisotropic linear theory itself, and its relation with possible non-analyticities of the effective shear moduli in 
the dilute limit. Bearing these considerations in mind, the present paper focuses on a new exact solution for the special, but 
representative, case of a two-dimensional (2D) array of voids embedded in an anisotropic linear elastic matrix, in the singular 
limit of infinite anisotropy. We emphasize that no exact solution of a similar type is available for general non-linear periodic 
media, which provides additional motivation for undertaking this study. 



B. Organization of the paper 



The constitutive laws, the prescribed loading conditions and the notations are defined in Sec. [Ill The word "loading" refers 
to the prescription of either overall conditions of homogeneous stress, or of homogeneous strain in the linear medium. Both are 
equivalent since they are related by the effective moduli. The matrix in the composite is compressible and possesses anisotropic 
properties in shear Two special limits of infinite anisotropy which lead to exact results are then considered (Sees. Hill and |IV|). 
For each of these limits, simple shear, pure shear, and equibiaxial loading situations are examined, and solutions are provided 
in each case (for some loading and anisotropy conditions however, only the limit of an incompressible matrix is considered; 
but this limit is the one relevant to non-linear homogenization of plastic porous media). The effective moduli and their dilute 
expansions, together with the mean and variance of the strain and stress fields are computed in each case. All these quantities 
are relevant to non-linear homogenization theories. A useful way of condensing the information contained in the solutions is 
by using field distributions which can in principle directly be used to compute the effective energy (Pellegrini, 2001) or, more 
trivially, the moments of the fields. Local extrema or saddle points in the field maps induce singularities in the distributions such 
as power or logarithmic divergences, or discontinuities, of the types observed in the density of states of a crystal (Van Hove 
1953; Abrikosov, Campuzano and Gofron, 1993). These singularities are examined in the context of our exact results. Cule and 
Torquato (1998) computed analytically the distributions of the electric field in a particular dielectric/conductor composite. They 
found singularities of the Van Hove type, but no extended singularities of the type put forward by Abrikosov et al. We extract 
below some singularities of this type, while singularities of yet a different type are encountered in Sec. |IVB"3] Obviously, in the 
periodic case, the field distributions are redundant with the exact solutions derived hereafter. However, for randomly disordered 
situations where exact solutions are not available, they remain the only way to collect useful informations on the fields. To ease 
their interpretation, a knowledge of their features in the periodic case is desirable, which provides a justification for the studies of 
Sees. nil B4l and |rV B 31 Moreover, we find that the infinite field variances obtained for some particular loadings are correlated to 
field singularities directly linked to localization patterns, which have counterparts as characteristic features in the distributions. 
A summary of our findings and a conclusion close the paper (Sec.lVll. 



II. PROBLEM FORMULATION 



A. Material constitutive law and microstructure 



We consider a periodic porous medium, with a square unit cell of size L made of a linear-elastic matrix (denoted by a phase 
index (3 = 1), containing a single disk-shaped void of radius a (phase index (3 = 2). Cartesian reference axes Ox, Oy are defined 
as in Fig.[T^, such that the void center lies at (x, y) = (0, 0). The unit cell is the square region [— L/2, L/2] x [— L/2, L/2]. The 
porosity is the surface concentration of voids / = 7r(a/L)^. The close-packing value where the voids touch is fc = 7r/4 ~ 0.78, 
for a = L/2. Hereafter, all lengths are rescaled by L, so that L = 1. The following notations are employed throughout the 
text: the upper-right quadrant (URQ) of the unit cell denotes the region {x,y) e [0, 1/2]'^, the lower-right quadrant (LRQ) 
stands for {x,y) G [0, 1/2] x [—1/2,0]. Similar abbreviations ULQ and LLQ stand for the corresponding left quadrants. 
Small deformations and plane strain conditions are assumed, so that the strain e derives from the two-dimensional (in-plane) 
displacement field u, and £y = if i or j = z. In-plane stress equilibrium requires diCTij = 0, i,j = 1,2 and we take the 
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FIG. 1: Left, periodic porous medium with unit cell and reference axes (a). Right, unit cell with 
(b) 0-type fibers; and (c) 45-type fibers. The black arrows depict eigenmodes of strain: simple 
shear (SS) in (b), and pure shear (PS) in (c). 



constitutive relation such that a^z = cTyz = 0, so that the problem is two-dimensional. The linear constitutive relation in the 
medium is cr(x) = L(x) : £(x). In the voids, L(x) = 0. In the anisotropic matrix phase, L(x) has components (Latin indices 
henceforth vary over 1 and 2): 

4Zi = ^'^-.^ + 2A Ef^% + 2fi EfJ^^i , (2) 

where k is the bulk compressibility modulus, where A, /i are in-plane anisotropic shear moduli, and where the operators JJ, E^^, 
E^^ are mutually orthogonal projectors of components (no summation over repeated indices here): 

A symmetric two-dimensional tensor a is expanded as a = I + ass ^ss + (^ps ^ps where I is the 2x2 identity matrix and 
where: 

ess = ( J J ) , ^PS = ( J ) , (4) 

the eigenvector sets of which are related by a 45-degree rotation. Corresponding strain modes are represented by arrows in Fig. 
[TJi and[T];. The equibiaxial component of a is a„i = [axx + o-yy)/'^- We call the second and third components the simple shear 
(SS) and pure shear (PS) components, with ass = ^xy and aps = {a^x ~ CLyy)/'^- Then J : a = a„jl, E^^ : a = ass ess and 
E^^ : a = aps eps- Accordingly, the constitutive relations decompose into: 

o-ss = 2Aess, crps = 2/ieps, CTm^ZKem. (5) 

Elastic tensor (|2]l is of a special type of orthotropy in the plane with Lim ~ L2222 7^ ^1122 7^ ^1212- When A > /i, this type 
of anisotropic response could correspond to a fiber-reinforced material with fibers aligned with the cell axes (see Fig.[TJ)). On 
the other hand, when > A, it would correspond to a material with reinforcing fibers aligned with the ±45° direction (see Fig. 
[1];). We introduce a shear anisotropy ratio, a, and normalized shear elastic moduli, m and t. 

a = A//i, m = /i/K, i = X/k. (6) 

The principal directions of the matrix being "aligned" with dense Unes of voids, a reinforcement of anisotropy-induced effects 
is expected. Besides, the loading modes considered below are aligned with the eigenstrains. These choices, motivated by future 
applications of this work to non-linear homogenization, focus on situations most relevant to plastic localization in porous media. 
Indeed, in the nonlinear theory, the loading determines the anisotropy of the background linear medium and always coincides 
with one of its eigendirections. Moreover, shear bands in random porous media quite generally link neighboring voids together 



B. Overall behavior 



In this work, we are concerned with the problem of determining the overall behavior of the periodic, porous material described 
in the previous subsection. This overall behavior is defined as the relation between the average stress (cr) and the average strain 
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(e) . Under the assumption of separation of length scales, the overall behavior of the porous material may be determined from 
the effective strain potential (see Torquato, 2002) 

M^=(l-/)min(,«W(£))^^^, (7) 

where denotes an average over the matrix phase, w'^'(£) = (l/2)£ : L'-^^ : e is the strain potential in the matrix and 

JC = {e\e — (1/2) [Vu + (Vu)-^] , u = (£)x + u*, u* periodic} is the set of kinematically admissible strain fields. The 
overall constitutive relation is then given by 

Because of linearity, we define the overall elasticity tensor L of the porous material via the relation (cr) = L : (£) . This tensor 
can be shown to take on the same form as dU with components Lij^ki defined by effective moduli A, Ji and H. 

In the analyses below, it will sometimes be more convenient to work with the effective stress potential U, such that the overall 
constitutive relation may be equivalently written 

(s)^l^, C/^(l-/)min(.«H)^^^, (9) 

where u*^^^(cr) = (1/2)(t : (L*^^-') ^ : cr is the stress potential in the matrix and S denotes the set of periodic stresses that 
are divergence-free in the unit cell, with prescribed average (cr), and traction-free on the boundaries of the pores. It should be 
noted that the case of an isotropic matrix (a — 1) has been addressed by McPhedran and Movchan (1994) in the more general 
framework of arbitrary contrast between the matrix and the square array of isotropic inclusions. 



C. Limits of infinite anisotropy and loading modes 

Hereafter only the strong anisotropy limits a — and a — +oo, amenable to an exact solution, are considered. Therefore: 
(i) When a — (i.e., X — or fi — oo) the medium is soft for SS loading, and resists PS loading; (ii) When a = oo (i.e., 
A = oo or /i = 0) the medium is soft for PS loading, and resists SS loading. Equibiaxial, SS and PS loading modes will be 
considered separately hereafter, so that only one of the three averaged components of the strain or stress is non-zero at a time. 
It is denoted by e (resp. a). Depending on the loading mode considered, strain and stress components and displacements enjoy 
various symmetry properties summarized in AppendixlAl 

For clarity, we detail the correspondence between the loading modes and the anisotropy properties of the material. Consider 
e.g. the limiting case a — X/ fi ^ 0, attained by two different types of materials: (i) A ^ and ^ > (Material 1); (ii) ii ^ oo 
and A < cxD (Material 2). We seek solutions with strain and stress finite almost everywhere (a.e.), i.e., except on points or on 
lines in the plane. In particular, aps — 2/ieps is finite a.e., so that eps = a.e. in Material 2. Likewise ess = fss/(2A) is finite 
a.e., so that ctss = a.e. in Material 1. Hence Material 1 is infinitely soft in the SS direction and can be finitely loaded in the PS 
direction, whereas Material 2 is infinitely rigid in the PS direction and can be finitely loaded only in the SS direction. Similar 
considerations apply for a — oo, leading to the following correspondence: 

a = 0: PS loading ^ A = 0, SS loading ^11 = 00, 
a = 00: SS loading ^ /i = 0, PS loading ^ A = 00. 

In equibiaxial loading, both types of materials should be considered for each of the limits a = 01 a = 00. 

As is well-known, the elasticity tensor (|2]i is positive definite when A, /i and k are all strictly positive, and we have existence 
and uniqueness of solutions. As mentioned above, however, the interest in this work is for the limiting cases where one of 
the eigenvalues of the elasticity tensor (|2]l (or, of its inverse, the compliance tensor) tend to zero. For these cases, care must 
be exercised when interpreting the solutions. As will be seen below, the relevant potential energy, or complementary energy 
functionals can still be shown to be strictly convex in the subspace of allowable strains, or stresses, and the standard theorems 
(see, for example. Proposition 1.2 of Ekeland and Temam, 1974) would still ensure existence and uniqueness of the solutions. 
However, appropriate components of the strain (or stress) field can develop discontinuities in these limiting cases. In this 
connection, it is also relevant to note that the governing equations lose ellipticity, leading to hyperbolic behavior (see below). 
This phenomenon is well-known in two-dimensional problems for ideally plastic materials (Kachanov, 1974), and has been 
exploited in a recent study of shape-memory polycrystals (Chenchiah and Bhattacharya, 2005). As was shown in this later 
reference, the hyperbolicity of the equations can be very helpful in interpreting the solutions obtained, and such connections will 
be made for the specific cases to be considered below. In any case, it should be kept in mind that the solutions of interest here are 
limiting cases of standard elasticity problems (with positive-definite elasticity tensors) for which solutions have been obtained 
by numerical methods (Willot et al., 2007). 
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D. Hyperbolicity conditions and characteristics 

Making use of the constitutive relations (|5]l, the stress equihbrium conditions may be written in terms of the displacement 
field as a closed system of two second-order partial differential equations (PDEs): 

+ a)^'"+(o f,ln)^>+^^ + ''-t'^{l J)a.9,u = 0. (10) 
Following Otto et al. (2003), this system is decoupled as: 

= 0, 0, Dcr-y = 0, i^x.y, j = x,y, (11) 

where the PDEs for the strain and stress fields follow immediately, and where D is the fourth-order differential operator: 

This PDE has the symbolic form (^{x, y) — x^ + y^ + 2rx^y^, with characteristic lines (e.g., Zachmanoglou and Thoe, 1986) 
defined by the equation ip{dx, dy) — 0, such that: 



dx ^ 
dy 



^J-r±{r^ - 1)1/2. (13) 



Taking A, /i, k > 0, r lies on the segment line r e [— l,cx)). At any interior point — 1 < r < oo, the quantity dx/dy is 
strictly complex; there are no characteristic lines in the real plane and the problem is elliptic. This corresponds to either (i) 
< A, /i < oo; (ii) < /i < oo, X = oo, k < oo; (iii) < A < cx), /i = cx), k < oo. At the end points r = — 1 and r — oo, the 
form (t){dx, dy) has two distinct real roots for dx/dy, each of them of multiplicity two. Thus the problem becomes hyperbolic 
with straight characteristics x/y = est: 

Case 1: r = oo <^ A = Oor/i = K = oo <^ x,y = cst, (14) 
Case 2: r = —1 ^ /i = or A = k = oo ^ x = ±y + est. (15) 

In this study, only hyperbolic problems are considered, i.e. the medium will be taken incompressible (k = oo) when A or 
11 = 00. All characteristics encountered are straight lines aligned with either one of the Cartesian axes (a = 0) or with one of the 
two diagonals of the unit cell (a = 00). As will be seen, the solutions found for the stress, strain or displacement fields actually 
verify simpler, first or second-order PDEs which can be used instead of (fTTT i. 



E. Averages, standard deviations and field distributions 

For any stress or strain component a(x) produced by a loading a, normalized phase averages and standard deviations in phase 
P are defined as: 

M^^Ha) EE {a)(p)/a, S''P\a) ^ ^{a^)^p^ - {a)f^^/a. (16) 

The distribution (i.e. the histogram) of a in the matrix M, Pa{t), is obtained by counting occurrences with the Dirac distribution 
5. All individual averages in the matrix result from the definition (g(a))(i) = j At g{t)Pa{t), where: 

Pa{t) = [l- f)-^ ( d^x5{a{-^)-t). (17) 

J M 



III. MATERIAL WITH ANISOTROPY RATIO a = 



A. Loading in pure sliear 

In this section, a = with A = 0. A PS stress loading a = (o-ps) is appHed. 
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FIG. 2: Structure of field patterns in the unit cell, in situations of infinite anisotropy. Figure (a): 
pattern for a = 0. Figure (b): a — oo (cf. Sec. lIVt . 



1. Stress fields 



Using (|5]l, A = implies ctss = a^y = 0. Then, from stress equilibrium, 

'^xx{x,y) ^ g{y), (jyy{x,y) ^ -g{x), 
where g is an unknown 1 -periodic function (using symmetry [ATJ. Hence: 

cTps.m = [g{y)±g{x)]/'2, 



(18) 



(19) 



where the plus (resp. minus) sign applies to ctps (resp. a^). This sign convention for PS loading, repeatedly employed hereafter, 
only holds for Sec. I/// A] The opposite convention will apply in Sees. \IIIB\ and \lirC\ 

We note that each stress component (Tu (no summation) obeys a simple first-order PDE, where x = est and y = est are 
characteristic lines, in agreement with (fT4l i. The stress vanishes in the void so that g{x) = for x G [—a, a]. The structure of 
the solution is most easily grasped referring to Fig.|2ja). It is organized in three types of zones: (A), (B), and (D), separated by 
"frontiers" marked by solid lines. The PS stress reduces to dps — g{x) /2 or g{y)/2 in zone (B), to crps = in zone (D), and to 
fps = [g{x) + g{y)]/2 in zone (A). It thus vanishes in the square of length 2a consisting in the union of (D) and of the void (V). 
This is more easily understood in terms of characteristics. Since the transverse stress component axy is zero, the void boundary 
conditions for the stress reduces to axx = o'yy = at any point of the void-matrix interface. Hence, each of these components 
is zero along one of the two families of characteristics (vertical or horizontal) passing through a void. At the intersection of 
characteristic lines passing through a void (i.e. region D), all stress components must be zero. Thus, as far as stress is concerned, 
the voids behave as square voids. In particular, the effective modulus (I34ab must be a natural function of a rather than of / (i.e. 
must not contain tt when expressed in terms of a). 

The unknown function g is obtained by minimizing the complementary elastic energy functional (|9]l, which is a strictly convex 
problem in the subspace of periodic stress fields with vanishing SS component (i.e., ctss = 0): 



-i<X,J/<i 

1/1 1\ /•1/2 
2 \/i K 



^crps{x,y) + ^cri{x,y) 



dx dy, 



g^{x) dx 



1 1 

fl K 



1/2 



g{x) dx 



Functionally extremizing this integral with respect to g under the constraint 

' ~ W=2 



(o-ps) 



1/2 

g{x) dx 



(20) 



provides g{x) = — 2a) for x E [—1/2, —a] U [a, 1/2]. Denoting by dixi,x2] (^) characteristic function of the interval 
[a;i,a;2],we introduce x(z) = 1 — 9[_g^g;^{z),the characteristic function of the domain [—1/2,— a]U [a, 1/2]. Then, crps and 
are completely determined by: 



""PS,: 



, X -x{y)±x{x) 
(x,y) = cr- 



(21) 



2(1 -2a) ■ 

Referring to Fig.|2ja), aps meets its highest (constant) value in zones (A), is of half this value in zones (B), and zero elsewhere. 
The discontinuous stress pattern obtained here as a solution for infinite anisotropy is of the type used by Drucker (1966) in his 
limit analysis of the ideally plastic porous medium. 
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2. Strain field 



From (|6]l and (1211 1. the strain components £ps and in the matrix read: 



/ N - x(y) - x{x) 



The relationship between cr and e = (eps) is now obtained. To compute e, we use u*. Expressions 
equation jAleU : 



1 a 

4/i 1 - 2a 



[(1 - m) x(a;) + (1 + m) x(y)]. 



(22) 

equivalently read [cf. 

(23) 



In PS loading, the admissibility (i.e., compatibility) conditions in (|7]i imply that dxU% = Exx — £ and dyU* — Syy + e. The 
displacement component u* (resp. u*) is odd and 1-periodic wrt. x (resp. y, see Eqs. lAll l. This requires u* to be tangent to the 
boundary of the unit cell: u*(±l/2, y) = u*{x, ±1/2) = 0. Hence, the previous Eqs. are integrated as: 



<(a;,2/) 



dx' [exx{x',y)-e], u*y{x,y) ^ I dy' [eyy{x,y') + e] . 

-1/2 J-l/2 



Furthermore, the periodicity condition m* (1/2, •) = yields: 



.1/2 

/ da; [exx{x,y) - e] = 0. 

J -1/2 



(24) 



(25) 



Choosing, e.g., y = 1/2 (in the matrix) and inserting (|23] ) in ( |25] ) provides the desired relation between a and e: 

_ 2^(1 - 2a) 
1 + (m - l)a 

Since m > and < a < 1/2, the denominator of ( |26] ) is always > 0. 

The still unknown Sxy = ess is computed from the admissibility (i.e., compatibility) conditions in Q, and from the expres- 



(26) 



sions of u* and u* in 



We obtain for (x, y) in the matrix: 



£xy{x,y) = - 



I dx' —Exxix' ,y) + I dy'—Syy{x,y') 
J-i/2 dy 7_i/2 ox 



(27) 



Inserting the solutions (l23l l. the integration is first carried out in the LLQ for integration paths in the matrix. The result is 
extended to the whole unit cell appealing to identities dAlfl g). For instance, in the URQ: 



ess.{x,y) 



1 + m 



4 1 + (m - l)a 



X - 5{y - a) ~ (y - ]^] 5{x - a) 



(28) 



where the Dirac distributions stem from the discontinuities in the displacements ( l24b due to the x functions in ( l23T l. Hence, in 
PS loading, ess is localized on the four "frontier" lines in the unit cell. Points (x, y) — (±a, 0) and (0, ±a) are singular hot 
spots. At each tangency point of a "frontier" line with the void boundary, there occurs a jump of egg along the line. E.g., in the 
vicinity of (a, 0), ( |28] | and the symmetry £xy{x, —y) = —Sxy{x, y) provide: 



£ss{x,y) ~ g^—T T^^ix-a) signy. 

8 1 + (m — l)a 



(29) 



From ( l24l i and identity u*{x, y) 
(9 is the Heaviside function): 

Ux{x,y) 



3. Displacement field and "hot spots". 
-w* (— X, —y), the displacement u* corresponding to the above strains reads in the URQ 



1 + (m - l)a 



(to — l)a + 2('^ + l)d{a — y) 



-(to — l)(a — x)9{a — x) 



(30) 
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FIG. 3: Structure of the deformed matrix and of u* (arrows) for an anisotropy ratio a = 0, in 
the particular case of equal bulk and shear moduli k. = /i. Loading is PS. Voids are in black. 
White (resp. grey) elliptic rings tag zones of extreme matter separation (resp. crushing). 



Moreover, u*{x,y) — —u*{y,x). Again referring to Fig. |2ja), u* is a linear function of x, y in (A), (B) and (D). Along 
the "frontiers" between (A) and (B) or between (B) and (D), its component normal to the "frontier" is continuous, whereas 
its tangential component is discontinuous. For instance, on the "frontier" y — a for a < x < 1/2, u* undergoes a jump 

[[<]]y = u*xix,a+) ~ u*{x,a-) given by: 

[[<]],(x,a) = -ui(l/2-a;), ^ (31) 

As expected, the displacement field can develop discontinuities along characteristic lines. This feature is reminiscent of rigid, 
ideally-plastic bodies in isotropic two-dimensional materials where discontinuities may only occur tangentially to slip lines 
(Kachanov, 1974, prop. 39.4). The structure of the full deformed configuration is schematized in Fig. [3) overall deformation 
occurs by block sliding. Fig.[3]is drawn in the particular case of equal bulk and shear moduli, k — ^ (m = 1). Then, u* vanishes 
strictly in zone (A) — but is finite there if k ^ ^ — and is aligned with the axes in (B). 
The expression of u in (D) clarifies the nature of the hot spots: 

u{x,y) =ui{signx,-signy), -a<x,y<a, (32) 

We stress that this incomplete expression, valid for ~a < x,y < a only, does not allow one to recover the singularity ( |29] l of 
£xy by taking a derivative. Since u is constant, save for orientation changes, in the four quadrants of zone (D), the hot spots are 
either points of extreme matter separation (at {x,y) — (0, ±a) for e > 0; white elliptic markings in Fig. [3]i or of matter crushing 
(at (x, y) ~ (0, ±a); dark markings). 

Fig. |3] also makes conspicuous four voided squares (in black) generated by the block-sliding pattern, at the intersections 
{x,y) = (±a,±a) of the "frontier" lines (the phenomenon is most remarkable for m = 1). The immediate vicinity of 
each of these points is divided into four regions where the displacement vector locally takes on distinct values. For instance, 
around {x,y) = (a, a), we have u*(a+, a+) = {uq, —uq), u* {a^ , a~) = (mq + Aw, — mq), u*(a^, a+) = (mq, — uq — Au), 
u* (a^ , ) = (uo + Au, —uq — Au), where 

Wo 2mi[(1 - m)/(l + m)]a(a - 1/2), Au = ui(l/2 - a). (33) 

The quantity Au represents the size of the voided square, to first order in e. Remark that in terms of a, it reads Au = 
ct(1 + m)/(8/i), a void-independent expression (the dimensioning factor is the cell size L — 1). 



4. Distributions, moments and ejfective shear modulus iJ.{f) 

In the PS case, the distributions Pe_ (t) or Per, (t) consist uniquely of a sum of Dirac-type components, plus unusual singular 
components at ±oo, due to the Dirac singularities in the fields. The latter cannot be accounted for by probability densities 
unless inconvenient limiting processes are employed. This may constitute a limitation of the use of distributions in a localizing 
regime. However, the first and second moments are readily computed. The effective shear modulus ]1 is read from (|26] l. With 

a = \/ f /n < 1/2, we have: 

- = 1 /r^"n =!-(! + ^)ifM'^' + (34a) 
^ 1 + (m — l)a 
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FIG. 4: Normalized shear effective moduli Jl/ and A/ A vs. void concentration /, for an in- 
compressible matrix with anisotropy ratio a = A//i = 0. Both moduli vanish at the close 
packing threshold / = /c = 7r/4 ~ 0.78. 



M^'\e,,)= M(^)(sps)^ V (34b) 

(1 - + (w- j[l + [m-l)a\ 



^ v/(l-2a)[(l + /)a-/] (1) , v/(l + /)a-/ 
^ (l-/)[l + (m-l)a] ' ^ ^""^- (l-/Vl-2a ' ^'"^^ 



[l + (m-l)a]VW' V(l-/)(l-2«)' 



(34d) 



5(i)(£ss) = oo, 5«('^ss)=0. (34e) 

The mean strain in the pore stems from A/ (2) (ep^) = [1 - (1 - /)M(i)(£ps)]//. Moreover, S*'!' (ctps) = {^/'fl)S^'^\eps)- The 
Dirac singularities are responsible for the transverse variance ^'■^•'(ess) being infinite. Also, remark that A/(^^(£ps) blows up 
as /^^/^ when / ^ 0. The curve pt(/)//i is displayed in Fig.|4]in the incompressible case m — 0. The power singularity at 
/ = with infinite negative slope, due to the /^/^ term in the dilute expansion of ( |34at , indicates that an infinitesimal void 
dramatically weakens the medium. This is reminiscent of the situation encountered in plasticity (see Introduction). In the linear 
material considered here, it is a direct consequence of hyperbolicity. Each of the two components axx and a-yy have constant 
values along one of the two families of characteristics. Due to the boundary conditions at the void-matrix interface, the parallel 
component crps is smaller by half in region (B) than it is in (A). Hence, each void lowers the stress field over large regions in 
the material, which are projections of the voids along the characteristic directions and involve infinitely long distances. The 
normalized standard deviations S''^-'(eps), S^-^\etn), S^-^\aps) and S'*^^-' (fJm) all behave as (//tt)^/"' as / ^ 0, which indicates 
that they grow more rapidly with / than in a linear isotropic medium. The effective modulus //(/) vanishes linearly with (/c — /) 
for / near fc ~ 7r/4, the void close packing fraction (see the comment in Sec. IIIIB 2b . along with the moments of the strain 
components in the matrix in the loading direction (A/^^^ (sps) — S^^^ (^ps) — 0). 



B. Loading in simple shear (incompressible case only) 

According to Sec. Ill CI SS loading goes along with a — and /i — oo. Matrix incompressibility renders the problem 
hyperbolic and is assumed for simplicity (k = oo, to = 0). 



1. Displacement, strain and stress fields 



Taking ^ = oo under finite stress implies e^x ~ £yy — 0. Incompressibility then requires Exx — £yy — 0. The only non- 
zero component is Exy, and a variational calculation is simpler to carry out wrt. the strain rather than to the stress. Condition 
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£xx = Syy = implies w* (x, y) = u%{y) and u*{x, y) = u*{x). By symmetry, see lA2cl u* (z) = u*{z) = G{z) — ez, where 
G is unknown. Introducing the derivative g = G , one obtains: 

e.,y{x,y)^\g{x)+g{y)\l2. (35) 

The energy of the unit cell ([7]l is expressed as an integral over the matrix phase M (represented by the unit square minus the 
void): 

=2A / d^x^^^^l / dxdy- / dxdyj[5(x)+5(2/)]2. (36) 

After using identity ( IA2fl ). we expand it into separate integrals over the intervals [0, a] and [a, 1/2], and we split g into indepen- 
dent functions g^ and g^ supported by these intervals. Subscripts A, B refer to the zones in Fig.|2a): 

9{z) ^ 9B{z)0[o^a]{z) + gA{z)9[a,i/2]{z)- (37) 
The strain energy ^ is functionally minimized wrt. gA, gs under the constraint: 

nl/2 



-xy/ 



gAiz)dz+ / 5B(z)dz 







(38) 



The system obtained determines gA as a constant. Introducing p{z) =\/ — z^, it provides an integral equation for gsiz)'- 

gA{z) = geia), zE [a, 1/2], (39a) 
2 I gB{y)dy ^ e + 2agB{a) + [2p{z)-l]gB{z), ze[Q,a]. (39b) 

Jp(z) 

In particular, for z = a, .gs (a) is expressed in terms of gsiz), z E [0, a[ : 



2/ dzgB{z)^e+{2a^l)gB{a). (40) 
Jo 

Taking z = in (|39b| l yields a relationship between g_B(0) and 93(0). 

gB{0) = [e + 2agBia)]/{l-2a). (41) 



Eqn. ( |39bt can be transformed into a second-order differential equation with no simple analytical solution. 1231 Instead we 
express the strain, the displacement and the stress in terms of 55, and we directly deduce from (I39bt some asymptotic results 
for the fields and for their distributions. We also compute numerically the whole solutions, the distributions, and the effective 
moduli by discretizing (I39bb into a linear system. The solution for gB{z) is represented v.s. < z/a < 1 for a E [0, 1/2] in Fig. 
|5] The inset displays .95(0) as a function of a. At a fixed, the highest values lie at z = 0, the strain being higher on the cartesian 
axes of the unit cell. Two regimes are encountered as a increases: the strain first develops in zone (B) and concentrates around 
the axes; then, for a > 0.4 (i.e. / > 0.5), (7b(0) blows up as (oc — a)~^/^ near Oc — 1/2 (see Sec. IIII B~3T l, while the strain 
localizes on the cartesian axes. In the limit, .95(2) oc S{z). The field enhancement between voids near close packing is boosted 
by anisotropy. The close packing void fraction = 7r/4 is, for periodic media, tantamount to the void mechanical percolation 
threshold in random porous media (e.g., Torquato 2002). When a Uc, gB{z) becomes discontinuous, with gB{o) = ^ 1 and 

55(0") = 0, see Fig. ID 

The solution in terms of gB is as follows. With (|39al l, g{z) in ( [37T i becomes completely determined in terms of gB{z), 
z E [0,a]. In the URQ: 



-^xy 



{X,y) = ^^Y^ [Ola,l/2]{x) + Ola,l/2]{y)] + ^[9B{x)0[o,a]{x) + gB{y)dlO,a]iy)'\- (^2) 



The integral equation (I39bl l admits a continuous solution. Then, gB and in turn Sxy, u and u* are continuous, too. The continuity 
and oddness wrt. x (resp. y) of u* (resp. u*) implies u* (0, y) = u*{x, 0) = 0. Hence, in the URQ: 

u*{x,y) ^ J dzg{z)-ey, u*y{x,y) = u*{y,x), {0 < x,y <l/2). (43) 
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z/a 



FIG. 5: Plots of the solution qb {z) of ( I39bb as a function of a/a, < a/a < 1, for values of the 
void radius (from bottom to top at a = 0) a = 0.025 (black), 0.124 (red), 0.249 (green), 0.275 
(blue), 0.373 (gray), 0.435 (violet), 0.485 (cyan), 0.498 (rose), 0.4995 (orange), 0.499985 
(dark green). Color online. Line and symbol thicknesses increase with a. 



These expressions are extended to the unit cell using (IA2aHA2c] i. From ( IA2db the transverse components (Jxx and ayy are odd in 
X and y, so that periodicity requires axxi^l'^^, y) = cryy{x, 1/2) = 0. With axy{x, y) — 2A exy{x, y) given by (|42] |. we deduce 
from stress equilibrium in the URQ of the matrix: 

.1/2 

(^xxix,y) = ayy{y,x) = / dx' dyaxy{x' ,y) = (A/2)(l - 2x)g'{y). (44) 

Then in this region: 

'JmM^, y) = (A/4)[(l - 2x)g'{y) ± (1 - 2y)g'{x% 

where the plus (resp. minus) sign applies to (Jm (resp. crps). Hence, because gsiz) = const, for z € [a, 1/2], the stresses cr,n and 
CTps vanish in zones (A) of Fig. |2 a). In zones (B), and crps are closely related, with: 

y) = apsix, y) = (A/4)(l - 2x)g's{y) for { ''<^<ll^ , (45a) 

y) = -apsix, y) = (A/4)(l - 2y)g'B{x) for { ^"4^/^;, . (45b) 

Taking the derivative of ( I39bl i with respect to z, we deduce 

p{z)[l ~ 2p{z)\ 

so that ,9^(0) = 0. Hence besides zones (A), regions of weak and dps lie along the axes a; = or y = 0, and along the axes 
a; = ±1/2, 2/ = ±1/2. 

Near the "frontiers" in Fig.|2ja), and crps blow up. Indeed, let us focus on the "frontier" y — a. From (l46T l. we deduce for 

g'siy) ~ -V2^[gB{a) + .gs(0)](a - y)-'/^ (47) 

and the latter stresses behave as d^^/^, where d is the distance to the "frontier" line. The variances involving an integral with 
integrand cx a^, this implies: 

S^^^Haps) = S^^\a^) = ^. (48) 
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These divergences are the only ones encountered in the means and variances. The mean stress readily derives from Eqs. ( [39pOl l 
and ( l42l l: 



1/2 r.l/2 r-a r.p(x) 

dx dy — dx dy 



Xs[l + gB{a)/e], (49) 



The effective shear modulus is therefore A = (A/2) [1 + (a)/£] . (50) 

2. Moments and effective shear modulus A(/) in the dilute limit 
Series expansions are now carried out for a ^ 1. Using the Taylor expansion: 

oo 

gsiaz) ^e^qn{z)a'^ (0 < z < 1), (51) 

we solve Eqn. (|39b| i. Note that a being a parameter, we should have denoted gsiz) by a), so that in the above equation, 
gsiaz) would stand for gsiaz; a), expanded in powers of a. The following recursion is obtained: 

qn+i{z) 



= gn(l) + Vl - z'^Qniz) - / ci/g„(y), qo{z) = l. (52) 

Numerically, we find that the series ( BTT l converges for a < Ra — 0.4. In the convergence region, this solution matches the 
numerical one obtained in the previous section. We carry out analytically the recursion up to fourth order to estimate gsiaz) /e. 
The integrations are readily performed with a symbolic calculator. Using the result with z = 1 provides: 

gs(a)/e = 1 - 27ra2 - + 2{Tr^ - Btt - 8)0^ + 0{a^). (53) 
With expansion (|53] |. we find in the dilute limit: 

^l-/-fS?+(l-^^)/'+o(/-). 

09 f3/2 09 fl/2 

5W(eps)=0, ^(i)(aps) = 00, ^^^^^^O, 5(i)K)=oo. (54d) 

A plot of A(/)/A obtained from the full numerical solution of gsiz) is displayed in Fig.|4] Compared to 'jl{f)/fi, the less 
singular character of the solution goes along with a harder material at small porosities. But A falls down to zero much faster than 
Jl for higher porosities, excepted near the close-packing threshold. 



3. Scaling in the close packing limit 

Day et al. (1992) examined the linear elastic behavior of an isotropic material containing a honeycomb lattice of voids: 
the effective compressibility and shear moduli (the overall medium is isotropic) vanish as (/c — /)^/^ near the void packing 
threshold. A similar behavior is found here in shear for the anisotropic material. 

With the numerical solution of ( |39b| ) near /c, we observe that the function gs (z) obeys a scaling of the form 



(55) 
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10 20 30 

z/(l/2-a) 

FIG. 6: Data collapse illustrating the scaling properties of the function (z), for the 3 different 
values of the radius a indicated in the legend, near the close packing value — 1/2. Inset: 
same plot in log.-log. scale; the curves break down when gB{az) become negative. 



where the master curve g{z) ^ const, for z ^ 1 and g{z) ^ for z ^ 1, where r is an exponent close to 2. Fig. |6]provides 
an illustration for three different values of a near Oc = 1/2. 

Two consequences are drawn from ( fSSl ): first, the scaling shows that the strain is concentrated in a band of width ^ ~ 
a(/c — jy^^, vanishing as the close-packing threshold is approached. Next, relation ( |55T l used at z = together with (HTb and 
(|50] | provides near f^. 

~>^/>^-{fc~ f Y'\ (56a) 

as for an isotropic material. The means are readily computed. Moreover, the variances are obtained with the help of the surface 
integral in the matrix J^^ da; Ay e^g {x,y) = e'^ [1 + (a)/e] /2. The following behaviors near the close-packing threshold ensue: 

m(i)(£ss) ~ (/c - ff\S^^\e^^) ^ (/e - ff'^ ., S^'Hass) - [fc - IT'" , 

^«(eps) = ^(^Hem) = 0, 5«(aps)-5«(a„) = oo. (56b) 

4. Distributions and Van Hove singularities 

The distributions of the stress fields in the matrix for SS loading, at void concentration / = 0.1, are displayed in Figs.|2]and 
[8]as thick solid lines. In these figures, straight vertical lines represent Dirac components. The distribution P^ss {t) of the stress 
"parallel" to the applied loading (Fig. [ill comprises one Dirac component, proportional to [(1 — 2a)^/ {1 — f )]5{t — gB{a)), which 
represents the contribution of zones (A) of constant (Tss. and a divergence followed by a discontinuity of infinite amplitude. The 
latter comes from the (B) zones. The right "foot" (magnified in the inset) is produced by zones (D) of maximal stress. The 
distribution has finite support since ess is bounded. The distributions of the transverse stresses cps and of are displayed in 
Fig. [8] They are even (these fields average to 0). A Dirac contribution at cr = (that of zone (A)) should be present in all plots 
but is omitted for clarity. The main components blow up at u = and slowly decrease at cr ^ oo (the contribution of zones (B), 
essentially). The distribution Pa„, differs only slightly from Paf^ by the contribution of zones (D), the difference being small at 
/ = 0.1, for the smallest values of the stress. 

The above singular features are understood as follows. Tails of the distributions at high values of the fields are determined 
by the square-root singularities responsible for the infinite variances ( |48T l. Definition (fTTT i of the distribution generates Van 
Hove singularities (VHS), in the vicinity of values t — Iq for which there exists extremal or saddle integration points Xj such 
that a(xi) = to and (9xa(xi) = 0. In 2D, depending on the eigenvalues of the matrix d^^a{xi) being of same or of opposite 
signs, the singularity in P is either a discontinuity or a logarithmic divergence (Van Hove, 1953). "Extended" Van Hove 
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FIG. 7: Distributions of the stress components parallel to the applied loading, at void concen- 
tration / = 0.1, in simple shear loading ((ass) = 1) with a = (Pass' solid), and in pure shear 
loading ((cps) = 1) with a = oo (Paps, dashed). Straight vertical lines near t = 0.8 indicate 
Dirac components. Inset: magnification of the "foot" in P^ss (t)- 




FIG. 8: Distributions of the stress components transverse to the applied loading, and of am, at 
void concentration / — 0.1, in simple shear loading ((ass) ~ 1) with a = (Po-ps (t) ~ Pa^ (t) 
for [t| < 1.5), and in pure shear loading ((crps) = 1) with a = oo (Passi^) — -Po-„(t)). Dirac 
components at t = should also be present in each of the four sets, due to the void and to zone 
A, but are omitted for clarity. 



singularities (EVHS) (Abrikosov, Campuzano and Gofron, 1993) arise when in addition one or more eigenvalues of d^^a{xi) 
vanish. Consider first the strain distribution. Using (|39b| i, we obtain 

9b{z)^9b{0)-^±^^z^ + O{z*), (57) 
a(l — 2a)'' 

so that in zone (B), for a<a::<l/2, 0<y<a, 

£ + ga(a) e + geia) 2 , ,.ox 
2(T3^-2a(l-2a)2^ +■■" ^^^^ 
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of the generic type h{x, y) = ho — hiy^, hi > 0. This generates an EVHS near t = Hq = {e + (a))/[2(l — 2a)], of the form 

Ph{t) ^ S{ho ~ hiy^ -t) = ^^^^{ho - t)-'^' 9iho - t), (59) 

where is the size of the integration domain on x values. Only values near y = contribute and the range of integration over 
y need not be completely specified. The four quadrants account for L^^ = 2(1 — 2a), so that 



P^^^it) c ^^MLj£^\t - hor^'9{ho ~ t). (60) 

Hence the divergence of Pa^s^ of square-root nature, and the drop of infinite amplitude in Fig.|7]are captured by this approxima- 
tion, since 

^.ssW=^e.„(i/(2A))/(2A). (61) 

Consider next the field distributions Pcrps [t) and P^^ (t) of Fig. [8] and examine the contributions near t ~ Q. According to (l45T l 
and to the remark following, weak field values lie on the axes x — 0, y — 0, x — ±1/2, y = ±1/2, where the four points 
(x, y) — (±1/2, 0) and (x, y) — (0, ±1/2) contribute most. Focus for instance on [x, y) = (1/2, 0). From (l45T l. ( |46] | and ( |49] |, 
and for y > : 

_ {l/2~x)y 

(Tm{x,y) ^ crps(x,y) ~ a— (62) 

a(l — Za)^ 



This is a VHS of type h{x, y) = hix^ - h2y^ (rotated 45°) with hi, /i2 > so that Ph{t) log \t\ / {2V y/hTh^) for t - 0, 

where V is the integration surface. Using ( |62] i, the contributions of the above four points gather into: 



' ' 2\t\ 



2(1 - 2a)^^ 



t 0. (63) 



The slowly-decaying tails at high t, of Paps{t) and Pcr„X'^) in Fig. [8] are due to singularity ( |47] i. Again focusing on the domain 
a < x < 1/2, 0<?/<ain zone (B), and expanding near y — a,we obtain from (|45] l , ( |46] l and ( |49] l : 

, A X / / N W^/2a 1 — 2x 

aps(x, y) ^ -(1 - 2x)g^{y) ^ -^^^-^^^==, y < a. 

Of type h{x, y) ~ — /io(l ~ 2a;) / y/a — y, Hq > 0, this expression contributes for 

lfa./"a.,M.„,-<)^|li^.H,.(l<,-li^), 

where the lower integration bound on y provides the restriction on |f | which delimits the domain of validity of the approximation. 
Adding contributions from negative values, and multiplying by 4 due to square symmetry leads to: 

This is only the contribution of (B). The tail should also include a contribution from (D) at high stresses, also decaying as \t\^^, 
but harder to handle due to sign changes at the points (±a, ±a). Other singularities exist which have not been fully investigated: 
values of around x,y = ±a/-\/2 produce a singularity in its histogram (extra wings starting at — ±1.7 with vertical 
slope, see Fig.[8]l; moreover, for porosities / higher than 0.6, two inflection points show up in 53(2). They give rise to two local 
extrema of and cps along the cartesian axes zone (B), also producing extra wings of the above type (not shown). 



C. Equibiaxial loading and limit A ^ 

We close the study of the case a = with equibiaxial loading e = (em) 7^ 0, (eps) — i^ss) — 0- The anisotropic limit a = 
follows either from letting A ^ 0, or /i ^ 00, so that two different solutions may a priori arise. Only the case A ^ 0, with 
finite K, is detailed in this section. Case /i ^ cx) is briefly addressed in Sec. 1111 Dl A method similar to that of Sec. 1111 Al for PS 
loading is used. The differences relative to the PS case can be related to the symmetries at play. 
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1. Stress and strain fields 



The limit A = under finite macroscopic mean stress requires that Uxy = 0. Equilibrium and the equivalence of the Ox and 
Oy axes imply that: 



(66) 



where g{z) is defined on the interval [—1/2, 1/2]. Similar steps apply as in Sec. IlII Al Energy minimization under the same 
macroscopic stress constraint (l20l i — with a = {a^) replacing (crps), provides (Jss{x, y) = 0, and 



o-m,ps(a;,y) = cr 



Xiv) ± X{x) 
2(1 -2a) ' 



(67) 



where the plus (resp. minus) sign applies to (Tm (resp. dps). Associated strains follow from Q. Periodicity conditions on u* 
entail: 



_ 2/i(l - 2a) 

= 

m + (1 — m)a 



(68) 



Eqn. (|67] | reproduces the PS result ( 1211 1. with the indices PS and m, and with x and y, exchanged. Furthermore, the PS and equib- 
iaxial loading modes exchange the parts played by k and /i, which corresponds to interchanging m ^ 1/to. This circumstance 
almost allows for a one-to-one mapping of the results of this section onto those of Sec. 1111 Al For instance in the URQ: 



Ess = U2 



1 



X - 7;]5{y-a)+\y- -]5{x-a) 



1 



1 + m 



U2 = — 



4 m + (1 — m)a 



(69) 



Eqn. ( |69l ) differs from ( |28] l by a substitution m ^ 1/m, and by a minus sign due to the difference between ( |Alg[ ) and ( |A3g[ ), or 
between dAlct and ( |A3c| ). 



2. Displacement field and singularities 



The displacement u* in the URQ is such that Ux{x, y) is there of the form dSOl l. with m replaced by 1/m, but with now 
Uyix, y) = u%{y, x). Hence u* is again piecewise linear, and tangentially discontinuous at the "frontiers" between (A), (B) and 
(D). In particular, the jump [[u^]] j, (a;, a) along the y = a "frontier" in the URQ is given by 0T1 ) with m replaced by 1 /m. With 
TO = 1, the overall pattern however differs from that of Fig. [3j by the fact that for e < (resp. > 0), the displacements u* in all 
parts of (B) are directed towards the void (resp. outwards), the medium being subjected to compression (resp. extension). 

Moreover, owing to symmetry u*{x, y) — u*{y, x), the PS mechanism responsible for the formation of the square voids at 
{x, y) = (±a, ±a) changes into one which induces a compressive singularity near these locations, independently of whether the 
equibiaxial loading mode is compressive or extensive, over a region of size Au = e{m + l)(a/2 — 1/4)/ [to + (1 — m)a\. This 
is readily seen from the following expression valid in the vicinity of (a;, y) — {a, a): 

u{x, y) = U2 [{a + 1/2) (1, 1) + (a - 1/2) (sign(y - a), sign(x - a))] . (70) 

The displacement in (D) now reads u(a;, y) — M2(sign x, signy), so that the four "hot spots" (x, y) = (0, ±a), (±a, 0) now all 
undergo a similar local singularity, either of a compressive, or of an extensive nature. 



3. Moments and effective compressibility modulus n 



The effective compressibility modulus and the moments are provided by expressions (I34ab -( l34el) . with /i and p replaced by k 
and K, with to replaced by 1/to, and with the indices PS and to interchanged. For definiteness: 



K _ (1 - 2a) _ 1 
/i TO + (1 — m)a m 



1 - 



1 



TO 



,1/2 



Oif) 



(71) 



Thus, K possesses a dilute-limit correction ~ /^/^ at finite k, but blows up as / in the incompressible limit m = 0. It decays 
as ifc - f) near f^. 
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D. Equibiaxial loading and limit /i ^ oo 

Consider now the limit oo, with finite A. A finite stress loading is consistent with this limit only if £ps = 0. This situation 
is that of a = under SS loading, the solution being of like complexity.ln particular, for finite compressibility, it is governed by 
elliptic equations (see Sec. Ill Dl ) and is much more complicated than the one of Sec. 1111 Cl Accordingly, we restrict ourselves to 
the incompressible case k = oo. The solution then simplifies and coincides with that of Sec. 1111 Cl with m = 0, see Note 
that K is infinite, due to /i = cx3. 

IV. MATERIAL WITH ANISOTROPY RATIO a = oo 

This section is devoted to the case a — oo, where the material is soft along the diagonals. As Fig. |2j5 illustrates, the field 
patterns are now rotated by 45 degrees. Accordingly, the solutions obtained below could alternatively be derived in a frame 
where the void lattice is rotated. The corresponding rotation, 72.450, would exchange in the matrix E^^ and E^^, and A and /x. In 
the rotated frame, the following correspondence should then be used 

72.450 (void lattice) <^ {a <^ l/a, and PS loading ^ SS loading, } . (72) 

Characteristics for the strain, stress or displacement fields are aligned with the two diagonals of the unit cell (see Sec. IllDI ). 
Another difference with the case a = resides in the appearance of new regions in the matrix, denoted by (C) in the figure, 
where the (B) bands cross, with no holes in it. 

A. Loading in simple shear 

We impose fj, = 0, a = 00 and consider SS loading conditions. Due to (|5]), /i = implies ups — {(Jxx — '^yy)l'^ = 0, so 
that axx{,'^^ y) = o-yy{x, J/) = s{x , Tj), an unknown function. Stress equilibrium implies d^s /dx^ — d^s /dy^ = 0. Using the 
identities in AppendixlAl we deduce that s{x, y) = [g{x — y) — g{x + y)]/2, where g is an unknown, even, 1-periodic function. 
The stress axy = crss stems from integrating the stress equilibrium equations. The integration constant is zero, since cr = in 
the void. Eventually, with g{0) — 0, 

crss{x,y)^[g{x-y)+g{x + y)]/2, a,^{x,y) = [g{x ~ y) - g{x + y)]/2. (73) 

Due to the void one has, along the main diagonal, amix, x) — —g{2x)/2 = for x < a/V2. Hence g{z) = on [— \/2a, V2a]. 
In turn, stresses ( |73] ) vanish in the whole square made of zones (V) and (D) in Fig.|2jb). Due to 1 -periodicity and to g being 
even, the stress pattern possesses a mirror symmetry with respect to the dashed lines of Fig. |2b). Its fundamental unit cell is 
therefore the gray square delimited by them. Symmetry ( |72] | makes this case tantamount to the PS case of Sec. 1111 Al with a = 
and A = 0, with a cell size diminished by a factor 1 / ^/2 at constant pore radius a and, consequently, with a porosity twice 
bigger. Effective moduli and moments are then obtained from expressions (|34|) of Sec. 1111 Al with fi and ]1 replaced by A and A, 
with TO replaced by £, with the subscripts SS and PS exchanged, and with a replaced by ^/2a. Deformation in the non-voided 
crossing zones (C) of zero stress is understood as follows. The obtained solutions for the displacement and strain fields in the 
(D) zones are trivially continued to the voided zone (V). These admissible continuations are non-physical except on the void 
boundary. However, the above symmetry considerations make clear that, modulo an appropriate translation, these continuations 
also provide the physical deformation fields in (C). 

Mutatis mutandis, the obtained expressions for the fields and the moduU are vaUd only up to / — fi^^ — it/8, due to the 
rescaling of a. This concentration acts as a "mechanically driven percolation" threshold, at half the close-packing value. It 
corresponds to the configuration where the summits (D) in Fig.|2]3 come into contact with the cell boundaries. Then, (A) and 
(B) vanish, leaving only (C) and (D) for / up to fc = tt/A. Thus, for / > fc^\ ass vanishes in the matrix. The stress ess also 
vanishes, so that strain only takes place along the void boundaries. Consequently, A vanishes for / > fc^^ and the fields are 
such that; M^^^sss) = 0, Af(i)(eps) = 1/f, S^-^\ess) = 0. Fig.|9]displays A(/)/A vs. /. Except for the threshold, the curve 
is alike that of Jl{f) / ^ in Fig. ID 

In words, we just showed that special field configurations in linear anisotropic periodic media can lower the geometric "per- 
colation" threshold of the porous material (i.e. the close-packing threshold, for a periodic pore lattice) into one determined by 
effective "porous" zones (as far as they undergo vanishing stresses) at intersections of void-generated stress bands. This is un- 
usual: for instance in homogenization methods, percolation-like thresholds are usually thought (or found) to be independent of 
the constitutive law. 
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B. Loading in pure shear (incompressible case only) 

In the case a = oo with X = oo and ^ finite, under PS loading, the solution has the pattern of Fig.|2tb). It is of the type studied 
in Sec lIIIBI at least for / up to /i^' = ir/S. For / > fj;^^ the problem becomes more involved due to a complex geometry, and 
has not been investigated in full. Only the range / < /i^^ is discussed hereafter. 



1. Displacements, strain and stress fields 

The solution obeys £ss = e^y = 0, and the incompressibility constraint e^^ + Syy = 0. Expressed in terms of m*, w*, such 
that u* (0, y) = u*{x, 0) = as implied by (lAlal i and ( lAlbk these equations entail with the help of (lAlck 

u*^{x,y) = [G{x + y) + Gix~y)]/2, u;{x,y) ^ [-G{x + y) + G{x ~ y)]/2, (74) 

and eventually with g{z) ~ G'{z) an even and 1-periodic function: 

£Ps(a;, y) = [g{y - x) + g{x + y)] /2. (75) 

Introducing the 45° counterclockwise-rotated coordinates x' — {y + x)/\/2 and y' — {y — x)/^/2, such that — 1/(2\/2) < 
x' Ty' < 1/(2 a/2), functions and gA are introduced, such that 

g{V2x) = 5B(a;')6'[o,Q](a;') +gA(a;')6'[^_i/(2v^)](a;')- (76) 
The total elastic energy in the unit cell ^ is provided by the integral: 

'J A'xel^^2J 2/ dx'dy' ^ [ dx'Ay'\el^{x\y'). (Jl) 

I- [-1/(21/2), 1/(2V2)]2 -IV J 

One extra contribution of the void V (of no elastic energy) has been subtracted from the contribution of the gray square in Fig. 
|2|b), counted twice. Anticipating on the fact that on its interval of definition g^iz) = g^{a), as in Sec. lIIIBTl the energy (|7]i is 
functionally minimized under the constraint: 

= 2\/2^gB(y)dy+ (l-2V2a) <?B(a), (78) 

This average is computed on one gray square, with the strain field ( iTSl l continued inside the void. Setting p{z) = \/ d?- — z^, an 
integral equation results: 

V2 l'^l^{y) dy = [1 - \/2p(z)l g^{z) - g^{a), z € [0, a]. (79) 
Jo '- 

Its solution is obtained following Sec. lIIIBTl In particular, with z = and dTSl l: 

5b(0) = ^+(l + 2V2a)gB(«). (30) 

' 2(1 -V2a) 

The displacement is everywhere continuous. The stress crps is deduced from (iTST i and from crps = 2/ieps in the matrix, with 
tTps = in the voids. The transverse components a^, css are computed analogously to (l44l i. with integrations and differentiations 
carried out over x' and y'. Stress equilibrium provides: 

dx' CTps ^ dy'{a,^- ass ) , dy' ctps = d^' {(Ttji + ass) ■ (81) 

Due to symmetry, a^nix, x) = ass{x, a;) = and crm(a;, 1 — x) — ass{x^ 1 — a;) = 0, so that with < a;, y < 1 in the matrix 
(for convenience, a unit cell translated by a vector t — (1/2,1 /2) is considered): 

a^,S5[x',y') = [l/2) [ dz dy'aps{z,y') ± (1/2) T dz d^'a^six' , z), (82) 

The plus (resp. minus) sign applies to am (resp. ass)- Then, in the matrix: 

o-m,ss(a;, y)^ ^[{x + y- l)g' {y - x) ± {y - x)g' {x + y - 1)] (83) 

for < X,?/ < 1. For z in [0, a\/2], the derivative g' reduces to g'{z) = g'^{z/^/2) /\/2 where g'g is the derivative of 
gB- The stress and strain singularities are qualitatively the same as in the a = 0, SS loading case. In particular, from ( i79l) . 
gs^z) ~ gsio^) + 2\/a[.9-B(0) + 9b{o.)]\<^ ~ for z < a, so that dps blows up along the "frontiers" of (A), (B), (C) and (D) 
as where d is the distance to the "frontier". 
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FIG. 9: Normalized shear effective moduli and A/A vs. void concentration /, for an 
incompressible matrix with anisotropy ratio a — oo. The curve for Jl/fj, is computed up to 



/ = 7r/8 only. Modulus A vanishes at /, 
value. 



(1) 



0.39, half the close packing threshold 



2. Moments and effective modulus fJ,{ f) in the dilute limit 
Expanding gsiz) in powers of a, we obtain: 

9Bia)/e = 1 - 27ra2 - ^^^s ^ [(5/2)^2 _ _ g] ^4 _^ Q(^^5y 
Dilute expansions in / for the effective modulus and the moments ensue: 



(84) 



(85a) 



M(i)(£ps) = 1 - + 0(f), Af(^)(£ps) = 1 + 



32/1/2 
3V27r3/2 



'0(/), 



(85b) 



V37r3/4 



(85c) 



(85d) 



In the range < / < 7r/8, is computed from a numerical solution of ( |79] l. Contrarily to the SS case where A(/) vanishes 
at 7r/4, no "mechanically driven percolation" is found here, since the zones (C) where bands cross are not zones of vanishing 
stress. Regular close packing behavior must occur at /c = 7r/4. 



3. Distributions 



Field distributions are studied as in Sec. lIIIB"4l Quite similar features are observed, summarized hereafter, along with some 
differences. The density P^ps is non-symmetric, and supported by the interval [g-^{a) , gB{0)]. A rescaling similar to ( |6T| l, using jj, 
instead of A, provides Po-ps, plotted in Fig. |7] (dashed curve). Due to the gradient of eps(a;, y) vanishing along the lines y — ±x, 
where the field is maximum in (B), Peps(^) blows up as \t — io|~^^^ at = [93(0-) + 5b(0)]/2. A notable difference with the 
case a = lies in the existence of the extended "foot" on the right side of P^rps ■ It ends with a discontinuous VHS of finite 
amplitude, due to the maximum strain field value, max Eps = 5b(0), being reached at {x,y) — (±1/2, ±1/2), in (C). In the case 
a = 0, <T = in this zone. The discontinuity is computed by remarking that according to ( |75] |, ( |76] | and to a Taylor expansion 
of gsiz) at z = 0, the strain field eps at points {x,y) — (±1/2, ±1/2) is a quadratic function with negative eigenvalues. 
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Furthermore, a field h{x,y) ~ —{hix'^ + /i2j/^) with hi, /12 > 0, on a domain V, produces in Ph{t) a jump at the origin 

Densities P^^, P^ss stem from Eqn. ( [83] l. They blow up at the origin due to points of vanishing derivatives dx<Tm,ss, 9j^cr„i,ss- 
According to ( |83] ) this happens for (x, y) = (1/2, 1/2). In the vicinity of this point, the expansion of is quadratic of a regular 
type and leads to logarithmic divergent VHS of Pa^it) at the origin. On the other hand, crss(a;, y) is of an unusual cubic form: 

ass{x,y)oc {y ~ x){x + y - 1) {x - 1/2) . (86) 

This generates a singular contribution to Pc^^^ (t) at the origin, of a type not previously encountered: indeed, the distribution 
Ph{t) of a function h{x, y) — ho{y^ — x^)x behaves near t = as jil^^/"^. Finally, the tails of Pcr,^{t), P„^^ (t) when |t| ^ cx) 
decay unsurprisingly as \t\~^, due to the blowing-up of the concerned stress components along the "frontiers" between zones 
(A), (B), (C), (D) as an inverse square root of the distance. 



C. Equibiaxial loading 



We consider here the limit a = A//i = 00, achieved by taking /i 0, under equibiaxial loading. The limit A ^ cx3 is 
not dealt with for the same reason as in Sec. IIIIDI Moderate porosities f < fc are assumed. The method used for the case 
a = is applied to the void lattice rotated by 45 degrees, and undergoing symmetry ( |72] |. With — \f2a and with the rotated 
coordinates x' — [y + x)/ \f2 and y' ^ [y ^ x) /a/2, the stresses are aps {x,y) = 0, and 

^m,ss(a;, y) = a[x{y') ± x(a;')]/[2(l - 2a,)], (87) 

in the rotated gray domain of Fig.lJJ). Hence, with w = {e/2)/[£ + (1 — £)as], 

£m{x,y) = iTuixiy') + x{x')], ess{x, y) ^ zu[x{y') - x{x')]- (88) 

Likewise, the PS component of the strain field takes the Dirac-localized form for x', y' in the URQ of the rotated gray square: 



eps{x,y) = {l+£)w 



V2 



5{y' 



Then, the following effective compressibility modulus is obtained: 

1 - 2a„ 



A 



£+ (1 -£)a. 



(89) 



(90) 



"Mechanically driven percolation" is again observed here. The moments of the fields are read from expressions ( I34bl i-( l34el i. by 
carrying out the following substitutions: (i) replace a by a,; (ii) replace to by 1 /t, (iii) replace index PS by index m; (iv) replace 
index m by index SS; (v) replace index SS by index PS. 



V. DISCUSSION AND CONCLUSION 



Two essentially different types of solutions, in an infinitely anisotropic linear elastic medium containing a periodic distribution 
of voids, have been exhibited. The anisotropy directions of the matrix — a "hard" one and a "soft" one — being aligned with 
the lattice directions, the fields are arranged in patterns of bands. 

On the one hand, for shear loading along a "hard" direction, the following observations were made: ii) the "parallel" compo- 
nents of the stress and strain are discontinuous, piecewise constant, with a band width of one pore diameter; {ii) the deformation 
pattern has a "block sliding" structure; (iii) the shear strain orthogonal to the loading has infinite variance, being localized as 
Dirac distributions along the band "frontiers", which are tangent to the voids (with jumps at the tangency points); iiv) the "per- 
pendicular" component of the stress is zero. The presence of discontinuities in the tangential component of the displacement 
field is analogous to similar observation in ideal plasticity, including the Hencky plasticity model (Suquet, 1981). Moreover, 
in situations where the bands cross, the medium develops fictitious porous zones of zero parallel stress which lower the close 
packing threshold. Also, a first-order correction ~ /^/^ is induced in the effective shear modulus. The latter decays linearly 
with the void concentration near the close packing threshold. 

On the other hand, for loading along a "soft" direction, another type of solution is produced when the problem is hyperbolic, 
i.e. assumed to be incompressible. In this case, the following observations were made: (i) the strains are continuous everywhere; 
(ii) the leading order correction to the effective shear modulus in the loading direction is now cx /, i.e. less sensitive to increasing 
void concentration for small concentrations; (iii) the variances of the "parallel" stress and strain in the matrix are cx p/'^ and 
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are again less sensitive to increasing void concentration; (iv) the parallel strain enjoys scaling properties near the close-packing 
threshold f^, operating within a band of width ^ ^ a(/c — fY^^, and becomes fully localized only for / near f^. This situation is 
at best one of "very soft" void-induced strain localization. However, stresses accumulate in the "hard" directions (transverse and 
equibiaxial) and blow up along the "frontiers" of the field pattern as the inverse square root of the distance (the corresponding 
strains being zero). This particular stress localization configuration may be relevant to strain-locking materials, such as the shape- 
memory polycrystals that have been considered recently by Bhattacharya and Suquet (2005) and Chenchiah and Bhattacharya 
(2005). 

In conclusion, situations responsible for fractional exponents showing up in anisotropic linear theories, of relevance to non- 
linear homogenization methods, have been clarified. Comparisons with fully numerical results at moderate anisotropics, and for 
isotropic viscoplasticity, will be presented elsewhere (Willot et al., 2007). As a final remark, the fact that strong singularities are 
found (isolated points of matter overlapping or of matter separation) in one case is obviously of relevance for breakdown studies 
(e.g., in brittle materials), since these points may act as initiators. This may suggest that mechanically sounder solutions could 
be looked for in a large deformation framework. However, in view of our more modest aim to help better understand possible 
hallmarks of localization in non-Unear homogenization theories, the model considered here is adequate. 
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APPENDIX A: SYMMETRY PROPERTIES FOR THE STRAIN, STRESS AND DISPLACEMENT FIELDS 

Since the periodic anisotropic medium obeys the symmetries of the square, the symmetries of the displacement u are deter- 
mined by the applied loading. They are easily deduced from Figs. [TJb) in simple shear and 1(c) in pure shear, and from an 
obvious flow pattern in equibiaxial loading. Differentiations wrt. x and y provides those enjoyed by . The symmetry group 
of flie constitutive law carries them unchanged over (Ty . For v = u or u* and for a = e or cr: 



• In PS loading: 



• In SS loading: 



• In equibiaxial loading: 
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